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Scattering at Q 2 =0.92 and 1.76 GeV 2 
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We report a Virtual Compton Scattering study of the proton at low CM energies. We have 
determined the structure functions Pll —Prr/e and Plt, and the electric and magnetic Generalized 
Polarizabilities (GPs) o.e(Q 2 ) and (3m(Q 2 ) at momentum transfer Q 2 = 0.92 and 1.76 GeV 2 . The 
electric GP shows a strong fall-off with Q 2 and its global behavior does not follow a simple dipole 
form. The magnetic GP shows a rise and then a fall-off; this can be interpreted as the dominance 
of a long-distance diamagnetic pion cloud at low Q 2 , compensated at higher Q 2 by a paramagnetic 
contribution from ttN intermediate states. 

PACS numbers: 13.60.-r,13.60.Fz 



The electric and magnetic polarizabilities of the nu- 
cleoli describe its response to a static electromagnetic 
field. Contrary to atomic polarizabilities, which are of 
the size of the atomic volume Q , the proton electric po- 
larizability ole is much smaller than one cubic fm, the 
volume scale of a nucleon. Such a small polarizability is 
a natural indication of the intrinsic relativistic charac- 
ter of the nucleon, as illustrated in a harmonic oscillator 
model The smallness of the proton magnetic polar- 
izability Pm relative to reflects a strong cancellation 
of para- and dia-magnetism in the proton. 

In Virtual Compton Scattering (VCS) -f*p — * 7p the 
polarizabilities become dependent on the momentum, or 
the four-momentum transfer Q 2 of the virtual photon, 
as first introduced by Guichon et al. 4] . These General- 
ized Polarizabilities (GPs) can be seen as Fourier trans- 
forms of local polarization densities (electric, magnetic, 
and spin) |5j . Therefore they are a new probe of the nu- 
cleon dynamics, allowing e.g. to study the role of the pion 
cloud and quark core contributions to the nucleon GPs at 
various length scales. VCS can be accessed experimen- 
tally via exclusive photon electroproduction ep — ► epj. 
After the NE-18 experiment @ and the pioneering VCS 
experiment at M AMI , we performed the E93-050 
H(e, e'p)j experiment [8| at the Thomas Jefferson Na- 
tional Accelerator Facility (JLab). We report low-energy 
expansion (LEX) analyses of our data up to pion thresh- 
old, and Dispersion Relation (DR) analyses of our data 
extending into the A-resonance region. 

To lowest order in the fine structure constant a em , the 
unpolarized ep — > epj cross section at small q' is: 



old photon production. The structure functions: 



d b a EXP = d 5 a BH+Born + q'Wo + 0(q' 2 ) , 
*o = vi ■ {Pll - ^Ptt) + v 2 ■ Plt 



(1) 



where 0, v\ , v 2 are kinematical coefficients defined in [tj] , 
q' is the final photon energy in the jp CM frame, and 
e is the virtual photon polarization. ( f j 5 a BH + Born corre- 
sponds to the coherent sum of the Bethe-Heitler (BH) 
and the VCS Born amplitudes, and depends only on the 
elastic form factors G P E ,G P M of the proton. This is a par- 
ticular case of Low's low-energy theorem ^(j f° r thresh- 



Pll - -Ptt 
e 



_ 4M p 




G E (Q 2 )a E (Q 2 ) + [spin-flip GPs] 



Pli = - — \I^G E (Q 2 )0 M (Q 2 ) + [spin-flip GPs](2) 

contain five of the six independent GPs [Hill- These 
structure functions are defined at fixed q, the CM three- 
momentum of the VCS virtual photon. Equivalently, Q 2 
in Eqs.© is defined in the q' — > limit: Q 2 — 2M p ■ 

(^M 2 + q 2 -M p ). 

The apparatus, running conditions and analyses of the 
JLab experiment are detailed elsewhere 0, ITil flH Il6l 
0,0]. An electron beam of 4.030 GeV energy was di- 
rected onto a 15 cm liquid hydrogen target. The two 
Hall A Spectrometers were used to detect the scattered 
electron and the outgoing proton in coincidence, allow- 
ing the identification of the exclusive reaction ep — > ep7 
by the missing-mass technique. This experiment makes 
use of the full capabilities of the accelerator and the 
Hall A instrumentation [l9|: 100% duty cycle, high res- 
olution spectrometers, high luminosities. We summarize 
our kinematics in Table ||| Variables such as q' , or the 
CM polar and azimuthal angles 6 and <p of the outgo- 
ing photon w.r.t. q, are obtained by reconstructing the 
missing particle. The acceptance calculation is provided 
by a dedicated Monte-Carlo simulation [2(j including a 
model cross section, resolution effects and radiative cor- 
rections plj . A number of cuts are applied in event anal- 
ysis, especially to obtain a well-defined acceptance and 
to eliminate protons punching through the spectrometer 
entrance collimator. 

We performed LEX analyses of the data sets I-a and II. 
The photon electroproduction cross section is determined 
as a function of q' , 9 and ip at a fixed value of q (1.080 and 
1.625 GeV/c) and e (0.95 and 0.88, respectively). The 
effect of the GPs on the cross section is small, reaching 
at maximum 10-15% below pion threshold. The method 
to extract the structure functions is deduced from Eq.([Q), 
in which the (BH+Born) cross section is calculated u sing 
a recent parametrization of the proton form factors [22j. 



For each bin in (0, <p), we measure d 5 a EXP in several bins 
in q' , and extrapolate the quantity AM — (d 5 a EXP — 
d 5 a BH+Born )/((j)q') to q' = 0, yielding the value of f . 
In our data, AM docs not exhibit any significant q'- 
dependence, so the extrapolation to q' = is done in each 
bin in (0,ip) by averaging AM over q' . The resulting 
term is then fitted as a linear combination of two 
free parameters, which are the structure functions Pll ~ 
P TT /e and P LT (Fig.HJ. 
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Q 2 =0.92 GeV 2 
-slope = +1.77 ± 0.24 
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FIG. 1: A graphical representation of the LEX fit (straight 
line) for data sets I-a and II. Circles correspond to out-of- 
plane data, and the inner plot is a zoom on the lepton plane 
data (triangles), ^o^i and V2 are defined in the text. 



The systematic errors are calculated from four sources 
added quadratically: 1) ± 2 MeV uncertainty in beam 
energy, 2) ± 0.5 mrad uncertainty in horizontal angles, 
3) ± 2.3% uncertainty in overall absolute cross section 
normalization, and 4) ± 2% uncertainty due to possible 
cross section shape distortions. The value of the reduced 
X 2 of the fit (Fig. Q is one measure of the validity of 
the LEX in our kinematics. The LEX results for the 
structure functions are summarized in Table ILTI 

We performed DR analyses of the data sets I-a and II, 
and also I-b including data from ttN threshold through 
the A resonance. In the DR formalism of Pasquini et 
al. [2^, the VCS amplitude is determined by unitarity 
from the MAID 7 (*)p — ► Ntt multipoles |24| , plus asymp- 
totic terms Act, A/3 which are unconstrained phenomeno- 
logical contributions to the GPs cce(Q 2 ) and (3m(Q 2 )- 



TABLE I: Kinematics of ep — + ep 7 . We used data sets I- 
a and II for the LEX analyses and all data sets for the DR 
analyses. 



data set 


Q 2 -range (GeV 2 ) 


IF-range 


I-a 


[0.85, 1.15] 


mostly < -kN threshold 


I-b 


[0.85, 1.15] 


mostly A(1232) resonance 


II 


[1.60, 2.10] 


mostly < -kN threshold 



Aa, A/3 are parametrized as follows: 
Aa(Q 2 )=a E (Q 2 )-a E N (Q 2 ) 



a 



exp 



,tJVi 



(3) 



(1 + Q7A 2 ,) 2 

(same relation for A/3 with parameter A^) where a E N 
(f3J^) is the 7rAf dispersive contribution evaluated from 
MAID, a e g P (/3^ p ) is the experimental value at Q 2 = 
0, and the mass coefficients A a and A^ are free pa- 
rameters. Theoretically, the choice of the dipole form in 
Eq.© is not compulsory. More fundamentally, the DR 
model provides a rigorous treatment of the higher order 
terms in the VCS amplitude up to the Ntttt threshold, by 
including resonances in the 7rA^ channel. In the region of 
the A(1232) resonance, these higher order terms become 
dominant over the lowest order GPs given by the LEX. 



TABLE II: Compilation of the VCS structure functions. In 
all cases the first error is statistical, and the second one is the 
total systematic error. 





Q 2 


e 


Pll - Ptt/c Plt 


(GeV 2 ) 




(GeV" 2 ) (GeV" 2 ) 


Ref. 






Previous experiments 


m 







81.3 ±2.0 ±3.4 -5.4 ±1.3 ±1.9 


® 


0.33 


0.62 


23.7 ±2.2 ±4.3 -5.0 ±0.8 ±1.8 


Set 




This experiment, LEX Analyses 


I-a 


0.92 


0.95 


1.77 ±0.24 ± 0.70 -0.56 ±0.12 ±0.17 


II 


1.76 


0.88 


0.54 ±0.09 ±0.20 -0.04 ±0.05 ±0.06 


Set 






This experiment, DR Analyses 


I-a 


0.92 


0.95 


1.70 ±0.21 ± 0.89 -0.36 ±0.10 ±0.27 


I-b 


0.92 


0.95 


1.50 ±0.18 ± 0.19 -0.71 ±0.07 ±0.05 


II 


1.76 


0.88 


0.40 ±0.05 ± 0.16 -0.09 ±0.02 ±0.03 



The DR analysis consists in fitting the free parame- 
ters A a and Ap to our cross-section data. This yields the 
value of the GPs a E (Q 2 ) and Pm(Q 2 ) using Eq.@. This 
also yields the value of the structure functions of Eqs.(|3 
since the DR model predicts all the spin- flip GPs |23| . 
Our DR results are presented in Tables ITTl and ITTll The 
systematic uncertainties are calculated from the same 
sources as in the LEX analyses. The error bars differ 
from one data set to another, due to differences in phase 
space coverage and in sensitivity to both the physics and 
the sources of systematic errors. The reasonably good \ 2 
of the DR fits (1.3 to 1.5) indicates that the DR model al- 
lows a reliable extraction of GPs in our kinematics, both 
below and above pion threshold. 

Figure [21 shows our DR extraction of the GPs cxe{Q 2 ) 
and (3m{Q 2 ), together with the point at Q 2 = [2( and 
the points derived from LEX analyses. The latter are 
obtained by subtracting the spin-flip polarizability pre- 
dictions |23] to the structure functions of Eqs.©- This 
involves some model dependence, which is not presently 
taken into account in the error bars. 

The solid curves in Fig. [3 are the full DR calculations, 
split into their dispersive 7rA^ contributions (dashed) and 
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The data for cxe(Q 2 ) disagree strongly with the simple 
dipole ansatz for the contribution Aa. It should be noted 
that our DR analysis is basically insensitive to the par- 
ticular choice of form of Aa and A/3, since our fits are 
performed independently in two small ranges of Q 2 . Fi- 
nally we point out that the r\ N and nirN channels, which 
must contribute to Aa, have resonances (S'ii(1535) and 
Di3(1520), respectively) with transition form factors that 
do not follow a simple dipole Q 2 dependence [2^, [2(j . 



o 0.5 1 1.5 2 

Q 2 (GeV 2 ) 



0.5 1 1.5 2 
Q 2 (GeV 2 ) 



FIG. 2: Compilation of the data on electric (a) and magnetic 
(b) GPs. Data points are from Ref. 0] (a), the LEX analysis 
of MAMI (o), and the present LEX (o) and DR (•) analyses 
of JLab. Some JLab points are shifted in abscissa for better 
visibility. The inner error bar is statistical; the outer one is 
the total error (statistical plus systematic) . Fhe curves show 
calculations in the DR model (see text). 



the remaining asymptotic contributions of Eq. J3J) (dash- 
dotted) for A Q =0.70 GeV and A^O.63 GeV, as fitted 
to the JLab data set I-b. The irN contribution to the 
magnetic polarizability in Fig. 0-b is strongly paramag- 
netic, predominantly arising from the A(1232) resonance. 
In the DR formalism, this is cancelled by a strong dia- 
magnetic term A/? originating from cr-meson i-channel 
exchange. The interpretation of A/3 as the contribu- 
tion of a long-distance pion cloud is further supported 
by the fact that the fitted scale parameter A^=0.63 GeV 
is smaller than the elastic form factor dipole parameter 
A=0.84 GeV. The dotted curves in Fig.[5]result from the 
full DR calculation, evaluated with Aq,=1.79 GeV and 
A^=0.51 GeV, which reproduces the MAMI LEX data. 



In summary, we studied the process ep — > epj at JLab. 
With data below pion threshold we applied the LEX, and 
for data extending through the A resonance we applied 
the DR formalism to extract the Generalized Polarizabili- 
ties. The different analyses are consistent, and the results 
give new insight into the correlations between spatial and 
dynamical variables in the proton. Other experiments at 
low energy will measure the VCS structure functions at 
low Q 2 [23, |28| and separate the six GPs via double po- 
larization measurements '28, 29] . 



TABLE III: The dipole mass parameters A a and A^ obtained 
by fitting the three data sets independently, and the electric 
and magnetic GPs evaluated at Q 2 = 0.92 GeV 2 (data sets 
I-a, I-b) and 1.76 GeV 2 (data set II). The first and second 
errors are statistical and total systematic errors, respectively. 



data set 


A Q (GeV) 


Ap (GeV) 


I-a 


0.741 ±0.040 ±0.175 


0.788 ±0.041 ±0.114 


I-b 


0.702 ±0.035 ±0.037 


0.632 ±0.036 ±0.023 


II 


0.774 ±0.050 ±0.149 


0.698 ±0.042 ±0.077 


data set 


a E {Q 2 ) (HT 4 fm 3 ) 


13m (Q 2 ) (10~ 4 fm 3 ) 


I-a 


1.02 ±0.18 ±0.77 


0.13 ±0.15 ±0.42 


I-b 


0.85 ±0.15 ±0.16 


0.66 ±0.11 ±0.07 


II 


0.52 ±0.12 ±0.35 


0.10 ±0.07 ±0.12 
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